AD-A089  651  SPERRY  RESEARCH  CENTER  SUDBURY  MA 

TINE “DOMAIN  SONAR  TARSET  RESPONSE  MOOELIN6.IU) 
OCT  79  H  MIERASr  C  L  BENNETT 
UNCLASSIFIED  SRC -CR- 79-74 


N61331-78-C-0049 

NL 


T 

i 


i 

4. 


MICROCOPY  RESOLUTION  TEST  CHART 


fldOO  31IJ 


05 

00 


TIME-DOMAIN  SONAR  TARGET  RESPONSE  MODELING 

\/  /  0*  L.feennett  - 

Sperry  Research  Center,  Sudbury,  MA  01776  * 


FINAL  REPpRI  . . _  J 

Contrap*-N«j  N61331-78-C-^049 J\r" 

i  ) 


(r,  'VS± 


RC-CR 


\-ivnL/ 


' jj  OctotaerT979  j 


Vv 


'■ir 

.  .fJTE 

SEP  2  9  1980 


D 


Prepared  for 

NAVAL  COASTAL  SYSTEMS  LABORATORY 
PANAMA  CITY,  FLORIDA  32407 


'DISTRIBUTION  STATEMENT  r 


Approved  for  public  reli 
Distribution  Unlimited 


^SPER^Y 

■  RESEARCH  CENTER 


JT-  6 


LIST  OF  ILLUSTRATIONS 


Figure 


Page 


2-1  Scattering  geometry.  4 

2-2  Frequency  response  of  sphere:  (a)  hard,  (b)  aluminum,  and  (c)  brass.  10 

2-3  Frequency  response  of  aluminum  sphere  with  (a)  ci  =  9300,  1 1 

(b)  ci  =  6200,  and  <c)  c1  =  4130. 

2-4  Frequency  response  of  aluminum  sphere  with  (a)  C2  *  4650,  12 

(b)  C2  =  3100,  and  (c)  C2  =  2065. 

2-5  Frequency  response  of  aluminum  sphere  with  (a)  p  *  2.7  and  (b)  p  *  1.5.  13 

2-6  Frequency  response  of  aluminum  sphere,  fluid  approximation,  with  14 

(a)  C!  =  9300,  (b)  c ^  =  6200,  and  (c)  c j  =  4130. 

2-7  Smoothed  impulse  response  of  aluminum  sphere:  (a)  hard,  (b)  aluminum,  17 

and  (c)  brass. 

2-8  Smoothed  impulse  response  of  aluminum  sphere  with  (a)  c]  =  9300,  18 

(b)  c1  «  6200,  and  (c)  =  4130. 

2-9  Smoothed  impulse  response  of  aluminum  sphere  with  (a)  C2  =  4650,  19 

(b)  c2  =  3100,  and  (c)  c2  *  2065. 

2-10  Smoothed  impulse  response  of  aluminum  with  (a)  p  -  2.7  and  (b)  p  =  1.5.  20 

2-11  Smoothed  impulse  response  of  aluminum  sphere,  fluid  approximation,  21 

for  (a)  ci  =  9300,  (b)  ci  =  6200,  and  (c)  ci  =  4130. 

2-12  Time  reference  for  scattering  from  unit  sphere.  22 


2-13  Aluminum  variations  -  responses  vs.  C2t  with  ci  =  6200.  24 

2- 14  Aluminum  variations  -  responses  vs.  C2t  with  c-j  =  9300.  25 

3- 1  Plane  wave  refraction  at  plane  boundary.  27 

3-2  Impulse  response  of  hard  sphere  (leading  edge).  32 


LIST  OF  ILLUSTRATIONS  (Cont.) 


Figure 

Page 

3-10 

Aluminum  response  model  —  time  domain. 

53 

3-11 

Aluminum  response  model  —  frequency  domain. 

54 

3-12 

Brass  response  model  -  time  domain. 

55 

3-13 

Brass  response  model  -  frequency  domain. 

56 

4-1 

STIE  response  for  fluid  aluminum  right  circular  cylinder. 

64 

4-2 

STIE  response  for  fluid  brass  right  circular  cylinder. 

65 

4-3 

STIE  response  -  small  angle  deviation. 

66 

4-4 

Surface  pressure  for  elastic  sphere  -  classical  and  STIE. 

70 

4-5 

Status  of  elastic  STIE  calculation  -  far  field  response. 

72 

4-6 

Shell  scattering  problem. 

72 

4-7 

STIE  solution  —  scattering  from  spherical  membrane. 

76 

5-1 

Shell  responses  as  a  function  of  thickness  —  frequency  domain, 
thick  shell. 

83 

5-2 

Shell  responses  as  a  function  of  thickness  —  frequency  domain, 
thin  shell. 

84 

6-3 

Shell  responses  as  a  function  of  thickness  -  frequency  domain, 
very  thin  shell. 

85 

5-4 

Shell  responses  as  a  function  of  thickness  -  time  domain, 
an  =  4,  thick  shell. 

86 

5-5 

Shell  responses  as  a  function  of  thickness  —  time  domain, 
an  *  4,  thin  shell. 

87 

5-6 

Shell  responses  as  a  function  of  thickness  —  time  domain, 
an  *  4,  very  thin  shell. 

88 

5-7 

Responses  of  very  thin  shell  with  variations  in  elastic  parameters, 
h-  0.001. 

90 

5-8 

Responses  of  very  thin  shell  with  variations  in  p  and  h,  keeping 
p^h  constant. 

91 

5-9 

Comparison  of  responses  of  soft  sphere  and  soft  core  (gas-filled) 
spherical  shell  -  frequency  domain. 

91 

5-10 

Comparison  of  responses  of  soft  sphere  and  soft  core  (gas-filled) 
spherical  shell  -  time  domain. 

92 

6-11 

Frequency  and  time  domain  responses  of  heavy  core  (glycerin-filled) 
spherical  shell. 

92 

SECTION  I 


INTRODUCTION 

The  acoustic  response  of  simple  metallic  targets  submerged  in 
water  is  studied  under  this  contract.  The  aim  is  to  develop  an  interpreta¬ 
tion  of  the  process  which  would  permit  its  characterization  by  a  simple 
model.  This  work  is  a  continuation  of  a  previous  contract  £>^7  which  should 
be  referred  to  for  the  general  background  of  this  problem  and  for  greater 
detail  of  the  space-time  integral  equation  approach. 


The  classical  solution  for  acoustic  scattering  from  a  penetrable 
homogeneous  sphere  is  studied  first  in  Section  2.  It  is  seen  that  hard  or 
soft  target  approximations  are  not  valid  for  metallic  targets  (except  at 
very  low  frequencies) ,  but  that  resonances  are  present  which  are  character¬ 
istic  of  the  interior  composition.  Viewed  in  the  time  domain,  the  impulse 
response  of  a  target  consists  of  a  sequence  of  pulses  following  the  specular 
impulse.  The  arrival  times  of  these  pulses  can  be  related  directly  to  the 
elastic  constants  of  the  target  and  medium.  Moreover  they  can  be  inter¬ 
preted  simply  as  pulse  traversal  times  at  the\several  speeds  along  certain 
paths.  The  description  of  this  model  (the  ♦glory-wave*  model)  is  the  sub¬ 
ject  of  Section  3,  along  with  the  first  steps  toward  a  complete  quantitative 
theory.  A  semi -quantitative  computer  program  wqich  generates  time  and  fre¬ 
quency  domain  responses  is  given  in  an  appendix 

In  Section  4,  the  space-time  integral  equation  approach  is  des¬ 
cribed.  Whereas  classical  solutions  can  only  be  found  for  targets  conforming 
to  separable  coordinate  systems,  the  integral  equation  approach  is  applicable 
to  target^  of  arbitrary  shape.  An  exact  solution  is  presented  for  a  fluid 
right  circylar  cylinder.  The  space-time  integral  equation  formulation  is 
presented  ftsr  elastic  targets  (which  exhibit  interior  shear  waves  as  well 
as  compression  waves) ,  although  numerical  difficulties  have  so  far  prevented 
computational  results.  A  procedure  is  described  whereby  the  scattering 
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fcenter  responses  can  be  computed  for  use  in  the  simple  model  of  Section  3. 

In  Section  5,  the  results  of  previous  sections  are  extended  to 
the  hulled  target.  The  classical  solution  is  derived  for  the  thick  spherical 
shell.  Results  fan  the  shell  are  given  as  functions  of  elastic  parameters 
and  shell  thickness  A  including  the  very-thin  shell  limit.  It  is  seen  that 
the  thick  shell  response  is  much  more  compl  i  already  compli¬ 

cated  homogeneous  sphere,  but  that  in  the  limit  of  the  thin  shell  the 
response  becomes  a  stepped  sinusoid  with  exponential  decay.  Also  given,  in 
subsection  4.3,  is  thei  space-time  integral  equation  formulation  for  the  thin 
hulled  target  of  general  geometry. 


/ 

/ 

/ 
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SECTION  2 

ANALVSIS  OF  HOMOGENEOUS  SPHERE  RESPONSE 


The  solution  for  scattering  from  a  homogeneous  penetrable  sphere 
was  first  given  by  Anderson  in  1950  for  the  fluid  sphere  [23  and  by  Faran 
in  1951  for  the  elastic  sphere  [33.  These  were  classical  solutions;  found 
by  expanding  interior  and  exterior  fields  as  sums  of  eigenfunctions  and 
equating  these  at  the  boundary.  The  applicability  of  the  method  is  limited 
to  simple  shapes  which  form  constant  surfaces  in  separable  coordinate  sys¬ 
tems.  However,  the  advantage  is  that  the  solution  is  straightforward.  The 
only  numerical  difficulty  arises  in  the  evaluation  of  the  spherical  Bessel 
functions  at  large  order  and  argument.  In  this  section,  the  classical  solu¬ 
tion  for  the  sphere,  transformed  to  the  time  domain,  will  be  used  as  a 
starting  point  for  understanding  acoustic  scattering. 


2.1  CLASSICAL  SOLUTION  FOR  HOMOGENEOUS  SPHERE 

The  elastic  sphere  is  completely  characterized  by  its  radius,  a, 

and  three  elastic  constants:  density  p^,  compression  speed  c^,  and  shear 

speed  c  (see  Appendix  8.1) .  It  is  embedded  in  a  fluid  with  density  p 
2  i  -*■  J 

and  speed  of  sound  c^.  Let  a  monochromatic  plane  wave  p  (r,t)  be  incident 

upon  it  from  the  positive  z  direction  (Fig.  2-1) .  The  incident  plane  wave 

can  be  expanded,  for  r  2  a: 

00 

«  P  V  (2n  + 1)  (-i)n  P  (cos  6)  j  (kr)  e~1U)t  .  (2-1) 

o  n  n 

n«0 

The  time  dependence  e~1Wt  will  be  suppressed  in  the  equations  to  follow. 
The  scattered  pressure  is  given  by 

00 

pS(r)  *  P  V*  c  h(1)  (kr)  P  (cos  d)  .  (2-2) 

r  o  Lmd  n  n  n 

n*0 
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Here, 


P^tcos  8)  are  the  Legendre  polynomials; 

j  (kr) ,  n  (kr)  ,  (kr)  are  the  spherical  Bessel  functions  of 

n  n  n 

the  first,  second  and  third  kind; 

Pq  is  an  arbitrary  constant,  or  could  represent  Pq(u))  ,  the 
frequency  content  of  the  incident  pulse; 

k  =  2TJ/X  =  w/c^,  the  wave  number;  and 

c  is  the  coefficient  which  is  to  be  found, 
n 

s 

Note  that  p  satisfies  the  radiation  condition,  since  asymptotically 

h*1^  (kr)  e  -*■  (-i)n+1  (1/kr)  -►  0  as  r  -►  00  and  is  outward  traveling, 

n 

-► 

Following  Faran,  it  can  be  argued  that  the  displacement  u  inside 
the  target  can  be  written 

u  =  -Vty  +  V  x A  ,  (2-3) 

where  \p,  A  are  of  the  form 

00 

*  =  E  an  jn(V)  Pn(cOS  6)  (2’4) 

n=0 

00 

A  =  $  ^  bn  jn(k2r)  ^  Pn(cos  6)  ,  (2-5) 

n=0 


with  unknowns  an  and  b^.  The  latter  two  equation  are  the  results  of  sym¬ 
metry  arguments.  These  expressions  are  now  evaluated  at  the  boundary; 
r  =  a,  and  related  via  the  boundary  conditions,  such  as  to  permit  solution 

for  the  unknown  coefficients  a  ,  b  ,  and  c  . 

n  n  n 

The  boundary  conditions  are  continuity  of  normal  stress ,  normal 
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displacement  and  tangent  shear  stress.  Or,  in  order,  in  spherical  coordi¬ 
nates  at  r  =  a: 


->•  3u 
A(V  •  u)  +  2y 

oT 


(2-6) 


9u 


0  = 


_e 

9r 


u 


_e  1  _ r 

r  r  96 


9u 


0  = 


r  sin  6  3<J> 


9u 

_ £ 

3r 


r 


(2-7) 


(2-8) 


(The  last  of  these  is  trivially  satisfied  due  to  axial  symmetry.)  See 
Appendix  8.1  for  the  relationships  between  the  Lame  coefficients  A,  u  and 
other  elastic  constants  and  for  the  relationships  between  pressure  and  dis¬ 
placement.  The  equation  of  motion  in  the  solid  is  (see  [1]  for  discussion) 


P1  u  =  ( A  +  2p)  V(V  •  u)  -yVxVxu  ,  (2-9) 

is 

so  that  (V  •  u)  and  (V  x u)  inside  as  well  as  p  =  p  +  p  outside  satisfy 
their  respective  wave  equations  with  \>elocities  c^»  c^,  and  c^. 

The  equations  (2-6)  through  (2-8)  must  be  satisfied  for  each  mode 
n.  We  thus  obtain  the  following  set  of  three  equations  (in  the  same  order) 
for  each  value  of  n,  by  direct  substitution  of  (2-1)  through  (2-3)  in  the 
boundary  condition  equations : 


l" 

a 

+ 

m” 

b 

+  M" 

c  =  v1 

11 

n 

12 

n 

13 

n 

i” 

a 

+ 

m” 

b 

+  m" 

> 

n 

o 

21 

n 

22 

n 

23 

n 

|" 

a 

+ 

Mn_ 

b 

*  0 

31 

n 

32 

n 

(2-10) 


n 

Mu  = 

|a(n2  +  nj  -  ■ 

xl)  JnK) 

-  2“  *1  j;(xl) 

n 

M12  - 

(n2  +n) 

o(x2 

j;(*2)  -  : 

>n(x2)) 

= 

6  D  a2 

h(l> 

(X) 

13 

o 

n 

n 

M  ,  = 

XX,  j  1 

(x\ 

21 

1  n 

V  1/ 

n 

M22  = 

(n2  +  n) 

*  jn 

W 

Mn  = 

2 

-Da 

h(1) 

’  (X) 

23 

o 

n 

n 

M31  = 

x,  j*  (: 
1  Jn  \ 

‘J- 

j„(xl) 

"32  ‘  -  4/2)  i„(x2)  *  X2 

n  ^ 

M  =0 
33 

Vn  *  -  6  D  (2n  +  1)  (-i)n  a2  j  (x) 

1  o  n 

Vn  =  D  (2n+l)(-i)n  a2  j’(x) 

2  o  n 

vn  =  0 

3 

and  where 


x  =  ka  =  0)  a/c^  , 

X1  =  kia  =  u)a^ci  ' 
x2  =  k2a  =  wa/c2  . 

It  was  useful  in  the  above  to  note  that  the  Bessel  functions  jn* 
h  all  satisfy 


(2-11) 


n  ,  and 
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r'l  - 


f"  (x)  =  (n2+n-x2)  f  (x)  -  2x  f  (x)  (2-12) 

n  '  '  n  n 


The  following  constants  were  defined  in  conformance  with  [1] 

2 


e  = 


P3C3 

2 

P1C1 


2c‘ 


a  = 


(2-13) 


D 


o 


We  can  now  solve  for  the  coefficients  a  , 

n 

M  .  This  must  be  done  for  all  modes  n. 


b  ,  c  by  inverting  the  matrix 
n  n 


By  trial  and  error  it  is  found  that  the  number  of  modes  n  =  NMAX 
required  to  achieve  convergence  is  about  NMAX  =  2  +  8v^a  for  low  values  of 
ka  and  about  NMAX  =  ka + 15  for  ka  greater  than  5.  The  applications  con¬ 
sidered  in  the  following  pages  required  ka  up  to  40  and  NMAX  up  to  55. 

The  evaluatic  of  Bessel  functions  for  such  high  order  and  argument  is 
non-trivial.  In  these  calculations  Univac  subroutines  were  used,  which 
were  based  on  the  work  by  Goldstein  and  Thaler  [4].  These  utilize  forward 

recurrence  for  n  and  backward  recurrence  for  j  .  There  are  still  diffi- 
n  n 

culties  with  these  evaluations  which  occur  when  n  is  large  but  the  argument 
is  small.  Overflows  in  the  computation  of  n^  then  limit  the  calcula¬ 
tions  to  lower  values  of  n  and  x.  This  situation  arises  when  c^c.^  is 
very  large.  For  the  materials  considered  here  (brass  and  aluminum),  it 
was  found  that  the  calculations  could  be  made  for  ka  up  to  60.  Since  the 


-8- 


effect  of  parameter  changes  was  also  considered,  the  calculations  of  fre¬ 
quency  responses  were  generally  limited  to  ka  =  40.  Further,  the  solution 
was  found  at  about  200  or  500  equally  spaced  points  on  the  frequency  scale. 

The  normalized  far  scattered  field  can  be  found  from  a  simplifi¬ 
cation  of  (2-2) ,  namely 


r  j  (r) 

— -  =  P  h(w)  , 

a  o 


c  (-i) 


n+1 


H  (u), 0)  =  \  '  ~n'  J~‘  P  (cos  0)  . 

- k -  n 

n=0  K 


(2-14) 


This  is  obtained  from  the  asymptotic  form  of  which  is 


h(1) 

n 


(kr) 


(-i) 


n+1 


ikr 


kr 


as 


kr  •+  00 


(Note  that  dropping  the  factor  e  in  (2-14)  is  equivalent,  in  the  time 
domain,  to  referencing  the  far  field  time  origin  to  the  space  origin  or 
sphere  center.  That  is,  a  pulse  scattered  from  the  origin  arrives  in  the 
far  field  at  time  t  =  0.) 


The  magnitude  of  the  frequency  response,  | H (to)  j  ,  in  the  back- 
scatter  direction  (9=0)  is  plotted  in  Figures  2-2  through  2-6.  Some  of 
these  results  were  essentially  previously  obtained  by  Hickling  [5].  (when 
making  comparisons,  see  the  note  on  "Scaling”  Section  3.8  in  [1].)  Responses 
for  sound-hard,  aluminum,  and  brass  spheres  are  given  in  Figure  2.2.  In 
Figures  2-3  through  2-6  the  responses  of  aluminum  are  given  and  compared 
with  their  modifications  due  to  parameter  changes.  This  will  be  discussed 
further  in  the  next  section.  It  is  seen  that  a  dominant  feature  of  the 
elastic  responses  is  the  regular  occurrence  of  deep  nulls  which  is  absent 
in  the  hard  target  and  fluid  target  approximations.  The  spacing  of  the 
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nulls  is  affected  by  the  velocities  c^  and  c^,  whereas  the  density  affects 
mainly  the  amplitude  of  the  response.  The  fluid  response  differs  from 
the  hard  approximation  mainly  by  amplitude  at  lower  frequencies;  although 
at  higher  frequencies  a  complex  structure  develops. 

These  features  are  difficult  to  explain  in  the  frequency  domain. 
After  all,  these  pictures  contain  only  half  the  available  information, 
the  phase  of  H(w)  is  as  important  (or  more  so)  as  the  magnitude. 

It  is  to  be  noted,  that  the  aluminum  and  brass  sphere  responses 
have  been  well  verified  by  measurement  [1]. 


2.2  SMOOTHED  IMPULSE  RESPONSE  FOR  SPHERE 

It  is  instructive  to  study  the  sphere  response  in  the  time 
domain.  To  do  this,  a  band  and  time  limited  pulse  is  chosen,  the  "smoothed 
impulse" : 

a  “(at')2 

e(t')  =  —  e  .  (2-15) 

/a 

This  pulse  has  been  normalized  to  give  unit  time  integral  with  the  normalized 
time  t'  =  c  t^/a  and  has  a  frequency  domain  representation 

-(ka/2a  )2 

E (ka)  =  e  n  (2-16) 

(Refer  to  [6]  appendix  for  a  discussion  on  normalization.)  In  equation 

(2-14),  P  is  replaced  by  E(ka) :  the  product  is  then  transformed  to  the 
o 

time  domain,  resulting  in  the  far  field,  normalized,  "smoothed  impulse 

S 

response",  rQ p  (r,t')/a.  The  larger  an»  the  sharper  the  smoothed  impulse, 

but  the  higher  the  range  of  gj  for  which  the  response  must  be  computed.  A 

practical  value  here  (consistent  with  the  considerations  of  Section  2.1) 

is  a  ■  8,  which  requires  the  computation  of  H(ka)  up  to  ka  *  40. 
n 
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2.3  EFFECT  OF  TARGET  PARAMETERS  ON  SPHERE  RESPONSE 

The  time  domain  responses  corresponding  to  the  cases  already 
given  in  the  frequency  domain,  are  plotted  in  Figures  2-7  through  2-11. 

The  time  axes  are  in  units  of  c^t/a. 

Consider  the  hard  sphere  response  of  Figure  2-7a.  The  initial 
part  is  the  specular  return,  an  impulse  at  time  t  =  -2.  This  is  immediately 
followed  by  a  negative  step,  a  ramp,  etc.  decaying  to  zero.  This  is  due 
to  radiation  from  the  region  immediately  surrounding  the  specular  point 
(see  Section  3.2) .  Then,  centered  at  a  time  t  =  it,  is  the  creep  response. 
Note  that  the  time  of  occurrence  is  just  that  of  a  pulse  inpinging  on  the 
target  at  z  =  0,  traveling  at  the  speed  of  the  medium  c^  around  the  back¬ 
side  of  the  sphere  and  reradiating  tangentially  at  z  =  0. 

It  is  seen  immediately  that  the  hard  target  cannot  be  an  ade¬ 
quate  model  for  scattering  of  wideband  radiation  from  an  elastic  target: 
there  is  a  small  change  in  the  specular  response  amplitude  and  backswing. 

But  mainly,  there  have  appeared  large  pulses  at  times  0  through  2.  (One 
can  verify  readily  that  interference  between  these  pulses  and  the  specular 
impulse  causes  the  resonant  structure  of  the  frequency  response.)  Signi¬ 
ficantly,  these  pulses  are  also  absent  from  the  fluid  model  response  (Fig. 
2-11) .  The  time  of  occurrence  of  these  pulses  is  strongly  affected  by  the 
shear  velocity  c^  (Figure  2-9)  but  only  slightly  by  the  compression 
velocity  c^  (Figure  2-8) :  we  thus  identify  these  features  as  due  mainly 
to  internal  propagation  of  shear  waves.  Also  present,  and  particularly 
noticeable  in  the  fluid  model  are  smaller  pulses  with  arrival  times  propor¬ 
tional  to  compression  velocity.  These  features  can  be  associated  with 
inter  1  compression  waves. 

Aiding  this  discussion,  Figure  2-12  illustrates  the  time 
references  for  various  reflection  paths.  Suppose  the  distances  a^,  a^, 
a3  are  travelled  by  a  plane  wave  at  velocities  c^,  c^,  respectively 


(as  compared  with  the  specular  wave  for  which  the  distance  is  a^  =  -2) . 
Then  the  time  of  arrival  is 


In  Figure  2-13  responses  are  plotted  versus  time  scale 

C2t  =  c2'/cl  al  +  a2  +  C2^C3  &3'  ^°r  severa*  values  of  c2"  If  the  Pat^  were 
entirely  at  speed  then  the  features  would  remain  at  the  same  value  of 
c^t.  This  is  nearly  the  case:  while  c^  is  changed  by  a  factor  of  4,  the 
time  c^t  at  which  the  "shear  wave  features"  appear  is  nearly  the  same. 

The  slight  variation  of  the  arrival  time  in  units  of  c^t  associated  with 
these  features  is  then  due  to  the  fact  that  the  paths  of  traversal  are 
in  small  part  at  c1  and  c^.  The  vertical  offset  in  Figures  2-13  and  2-14 
is  proportional  to  .  Figure  2-14  is  the  same  as  2-13  but  with  a  higher 
value  for  c^.  Only  for  large  values  of  c 2  is  there  a  significant  difference 
between  the  curves.  These  results  are  sufficient  to  allow  estimating  path 
lengths  through  the  target  at  the  various  velocities  for  the  principal 
features  of  the  response.  Physical  considerations  (Snell's  law  —  Section 
3.1)  then  can  give  the  likely  path  associated  with  these  features.  In 
Section  3.3  they  will  be  identified  as  the  principal  modes  of  a  glory-wave 
effect. 
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c2/e20 


SECTION  3 


SCATTERING  MODEL 


It  is  apparent  from  both  frequency  and  time  domain  points  of 
view  that  the  scattering  process,  for  even  such  a  simple  target  as  a  sphere, 
is  inherently  complicated.  The  aim  of  this  section  is  to  identify  the 
essential  attributes  of  the  responses  in  order  to  obtain  a  reasonable 
approximation  with  a  relatively  pimple  model.  Such  a  model  is  possible 
based  on  the  glory-wave  effect.  This  theory,  developed  originally  for  the 
scattering  of  light  from  water  droplets,  undergoes  considerable  complica¬ 
tion  (due  to  the  existence  of  two  sound  speeds  in  the  target)  when  applied 
to  the  acoustic  case.  However,  the  qualitative  aspects  of  the  theory, 
particularly  when  viewed  in  the  time  domain,  are  helpful  in  providing 
useful  insight  into  the  scattering  process;  it  will  be  seen  that  the  im¬ 
pulse  response  consists  of  a  string  of  inpulses  which  arrive  after  trav¬ 
eling  through  the  target  along  the  readily  identifiable  glory  paths.  On 
the  other  hand,  quantitatively  (i.e.,  for  the  magnitude  of  the  impulse 

response)  the  theory  is  difficult  to  apply  (see  the  extensive  work  of 

•  • 

Uberall  etal.  [7].  It  is  suggested  that  a  useful  approach  is  a  hybrid 
approximate  method,  using  the  quantitative  considerations  of  the  next  two 
subsections  together  with  the  simple  glory-path  model. 


3.1  REFLECTION  AND  REFRACTION  AT  A  PLANE  ELASTIC  BOUNDARY 

The  following  is  an  exact  treatment  of  the  interaction  of  a 
plane  compressional  wave  with  a  plane  fluid-elastic  interface,  and  of  the 
reverse  problem. 

Let  the  displacement  in  a  fluid  medium  be  given  by  u  ,  and  in 
the  solid  by  u^  and  u^  corresponding  to  compressional  and  shear  waves 
respectively.  Primes  are  used  to  denote  the  part  of  the  wave  after  reflec¬ 
tion  or  refraction.  Consider  the  plane  interface  of  Figure  3-1.  Then  we 
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(b)  GENERAL  CASE 

7»-«’0 

FIG.  3-1  Plane  wave  refraction  at  plane  boundary. 


have 


incident  wave  =  u^  e 


jk3  •  r  -  joit 


(3-1) 


reflected  wave  u*  =  u^Q  e 


jk^  •  r  -  jo)t 


(3-2) 


refracted  wave  u'  =  u*  e 
1  10 


jk'  •  r  -  jwt 


(3-3) 


U2  U20  6 


jk^  •  r  -  juit 


(3-4) 


The  zero  subscript  refers  to  the  value  at  the  origin.  Assume  that  the 
wave  normals  are  in  the  x,z  plane  so  that: 


k3  =  (x  sin  93  -  z  cos  03)k3 

k^  =  (x  sin  03  +  £  cos  03  j  k^ 

k|  *  (x  sin  6^  -  z  cos  0^J  k^ 

k^  *  (x  sin  02  -  z  cos  02]  k^ 

where 

k.  =  u)/c.  . 

i  i 


(3-5) 


Snell's  laws  for  reflection  and  refraction  follow  directly  from  these 
definitions,  independent  of  the  nature  of  the  boundary  conditions  (see 
Jackson,  p.  217  [8]) :  whatever  boundary  conditions  prevail,  the  same  con¬ 
ditions  are  satisfied  at  the  point  (0,0)  as  at  every  other  point  on  the 
interface.  Hence  the  variation  in  phase  for  (3-1)  through  (3-4)  must  be 
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the  same  for  all  under  translation  on  the  interface.  That  is: 

sin  0^  =  sin  0 =  k^  sin  0^  =  k3  sin  0^  (3-6) 

(The  same  argument  for  phase  variation  with  time  can  be  used  to  establish 
that  the  frequencies,  w,  are  the  same  for  all  waves.)  Hence  Snell's  Law: 

sin  6^  =  sin  0^ 

c! 

sin  0  =  —  sin  0,  (3-7) 

1  C3  3 

C2 

sin  0  *  —  sin  0_ 

2  c3  3 

Note  that  for  solids  we  generally  have  >  c2  >  c3»  so  that  rays  are 

refracted  away  from  the  normal  after  entering  the  solid,  as  in  Figure  3-1. 

Also  note  that  for  some  real  incident  angle,  03cl  (the  critical  angle),  we 

will  have  sin  0,  =  c,/c„  sin  9,  ,  =  1.  There  is  no  transmission  at  speed 

lc  1  3  3cl  * 

c,  for  0„  2  0„  , .  Similarly,  at  0,  _  such  that  sin  0_  =  c  /c  sin  0„  _  = 

1  3  3cl  3c2  2c  23  3c2 

there  is  no  transmission  at  c 2  (nor  at  c^) ;  so  that  there  is  total  reflec¬ 
tion  for  0  i  6  At  angles  between  these  two  critical  angles,  there 

is  transmission  only  at  c^,  that  is  as  a  shear  wave.  The  boundary  condi¬ 
tions  can  be  'satisfied  through  the  existence  of  an  evanescent  wave  which 
travels  along  the  surface. 

We  will  now  derive  the  relative  magnitudes  of  these  waves, 
for  which  it  is  necessary  to  know  the  dynamics  or  boundary  conditions. 
These  are 

(1)  continuity  of  normal  stress, 

(2)  continuity  of  normal  displacement, 

(3)  continuity  of  tangent  shear  stress. 
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In  Cartesian  coordinates  these  are  written  [1]: 


x3[(v-;3)  ♦  (?•:•)]  -  a<v •  ;>  ♦  2u  -g  , 

n  •  [u3  +  ».]  =  n  •  u  , 

f9u  3u  *1 
+ 


where  u  =  +  u2  and  the  elastic  constants  are 


VP3  =  °3 


(X  +  2y)/p  =  c1 


M/P  =  c2  . 


(3-8) 


(3-9) 


Substituting ' equation  (3-1)  etc.  into  the  boundary  conditions,  we  obtain 
respectively 

Sk3  (U30  +  u30>  '  Xkl"io  +  2U[uio  kl COs2  91  *  U20  k2 cos  62  Si"  62]  ’ 

(3-10) 

(u30  ‘  u3o)  COS  ®3  *  U10  cos  61  *  U)o  Sin  62  • 

2  u30  k3  sin  03  cos  0^  -  *2  (cos2  02  -  sin2  ^  j  =  0  .  (3-12) 


Solving  simultaneously  we  obtain  the  transmitted  amplitude  (that  is  u^, 

u'  ,  and  u'  )  in  terms  of  the  incident  amplitude  u  n.  It  is  convenient 
20  30 

to  define  the  quantities 


30- 


=  ^|l  -  a  sin2  0^1  +  aQ  sin  0^  cos  02^ 


(3-13) 


cos  0. 


B  = 


(2  2 
cos  0„  -  sin  0_ 
2  2 


) 


(3-14) 


Q  = 


2  c^  sin  0,  cos  0, 

2  1 _ 1_ 

c, /cos2  @2  -  sin2  02j 


^(C 


(3-15) 


2  2 

where  a  =  2  c^/c^.  Then  the  transmission  coefficients  are  found  to  be 


R  = 


u* 

30  A  -  B 


33  U3Q  A  +  B 


(3-16) 


10 


3!  U3Q  A +  B  ' 


(3-17) 


R  =  _20  ,  JSL. 

32  u3Q  A  +  B  ’ 


(3-18) 


The  first  subscript  of  R  refers  to  the  incident  wave,  the  second  to  the 
reflected  or  refracted  wave. 


The  above  derivation  can  be  performed  more  generally,  as  indi¬ 
cated  in  Figure  3-2,  to  also  permit  waves  incident  from  inside  the  solid. 
The  results  are  as  follows.  For  an  incident  compression  wave  in  the 
solid,  u1Q,  the  refractions  and  reflections  are 


'll  u 


12  u 


uio  _ 

B  -  C  +  QD 

U10 

A  +  B 

u’ 

20 

11  ■  *S 

.  J2c_ 

U10 

A  +  B 

(3-19) 


(3-20) 
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FIG.  3-2  impulse  response  of  hard  sphere  (leading  edge). 


_  ^30  =  2BC 
^13  '  u1Q  A+B 


(3-21) 


Here  we  have  also  used  the  definitions 


£ifi  / 

p^s  \ 


COS2  “*  s^n  ^9  |  ' 


(3-22) 


P1C2  ( 

°  P,c,  v 


2  sin  02  cos  ©2 


(3-23) 


the  refractions  and  reflections  are 

For  an  incident  shear  wave ,  u2q  * 


10  _  2P 

R21  u2Q  A  +  B 


(3-24) 


(3-25) 


22  u„ 


_  ^30  =  _  2BP 
*23  "  u2Q  ”  A  +  B 


(3-26) 


Note  that  R13.  R23  represent  energy  transmit 
while  R._  and  R  represent  mode  conversions. 


transmitted  from  solid  to  fluid  medium. 


normal  incidonce.  \  -  «2  -  «3  -  ».  .xprossions 


plify  to 


_  a^i  .  pici :  p3c3 

*33  '  A  +  1  +  P3C3 


2P3C3 


*31  A  +  1  Pj^  +  P3C3 


(3-27) 


R  *  0 
32 


Also,  it  is  interesting  to  investigate  the  case  of  critical  incidence  for 
c1»  that  is  =  90°.  We  then  have  that  R33  =  1,  that's  total  reflection. 
As  the  incident  angle  increases  slightly  and  0^  becomes  imaginary,  there 
is  transmission  of  energy  as  a  shear  wave  at  c^,  while  on  the  surface 
there  exists  an  evanescent  compression  wave  at  c^. 

In  closing,  it  is  noted  that  the  fluid,  hard,  and  soft  cases  are 
included  in  the  analysis  by  appropriately  letting  the  densities  p^,  p 3 
go  to  zero  or  00 . 


3.2  TIME  DOMAIN  SPECULAR  RESPONSE  -  LOCAL  STIE 


In  the  physical  optics  limit,  the  impulse  response  of  a  target 
is  given  by  1/2tt  times  the  second  derivative  of  its  projected  area  func¬ 
tion.  This  was  given  by  Kennaugh  &  Cosgriff  in  1958  [9]  for  the  electro¬ 
magnetic  case  and  is  valid,  as  will  be  seen,  in  the  limit,  for  the  acous¬ 
tic  case  as  well.  For  example,  the  projected  area  of  a  unit  sphere  with 
specular  point  at  the  origin  (Figure  2-1)  for  an  incident  plane  wave 
traveling  in  the  z-direction  is 


(3-28) 


Taking  the  speed  of  sound  c  =  1,  we  have  that  z  =  t/2.  The  area  function 
is  then 


S  (t) 
P 


(3-29) 


and 


(3-30) 


is  the  physical  optics  impulse  response.  This  can  be  considered  to  be  the 
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first  two  terms  of  an  expansion  in  singularity  functions.  An  improvement 
to  this  approximation  can  be  obtained  by  including  the  effect  of  surface 
interactions . 

3.2.1  Sound  Hard  Convex  Specular  Point 

For  a  sound  hard  target  the  far  backscattered  pressure  is  given 
by 


ro  pS(t) 


*') 


3p(r',f) 

3t* 


dS' 


(3-31) 


where  t'  =  t  +  (r  •  r')  and  r  is  the  direction  of  observation.  The  pres- 
o  o 

sure  on  the  surface  can  be  shown  [1]  to  be 


P(r\f) 


2  p  (z’,f) 


(R 


i  -JL\ 

R  9t"  j 


p  (r  " 


t”  )  dS"  , 
(3-32) 


where  t”  =  t*  -  R,  R  =  r  -  r' ,  and  t  is  in  units  of  sound -meter s .  It  is 

easier  to  consider  the  step  response  first.  Let  the  incident  plane  wave 

A  ■+ 

be  a  step  traveling  in  the  positive  z-direction  (i.e.r  rQ  •  r'  =  -z') 


P  (z'.f) 


u(t’  -  z ' )  , 


or 


p1(z',t)  =  u ( t  -  2z ' )  .  (3-33) 

•f 

As  a  first  approximation  we  will  consider  the  surface  pressure,  p(r"  ,  t"), 
constant  over  the  illuminated  region  (and  zero  outside) .  Then  3p/9t"  =  0, 
except  at  the  boundary  of  the  illuminated  region  where 
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7,  (r"  ,  t">  =  p(r"  ,  t")  6(t"  -  z "  )  , 


(r" ,  t")  =p(r",t")  6(t‘-zn  -  R)  . 

We  want  to  substitute  this  in  the  integral  (3-32).  To  evaluate  the  inte¬ 
gral  we  use  the  fact  that  for  a  circular  patch  As  of  radius  Y 


•£/ 


...  LS  *  _  £ 


,  Y< 
dR=  2 


(3-34) 


where  <  is  the  average  curvature.  This  approximation  applies  to  smooth 

convex  surfaces.  This  was  derived  in  [1]  p.  16,  but  due  to  a  typographical 

error,  the  value  of  £  is  given  incorrectly  in  [1].  Here,  As  is  the  ellip- 
— ► 

tical  patch  around  r'  at  time  t'  which  is  illuminated  at  the  retarded 
times  t"  =  t'  -  R.  Some  reflection  reveals  that  the  approximate  radius 
of  AS  is  given 


Y  =  t  -  2z'  . 


For  the  second  term  in  (3-32)  involving  3p/9t"  we  need 


(3-35) 


*/ 


ds..  -(-  ’  -n— -  6(t'  -  z"  -  R)  =  J 


I  _ 

/  R  6(t*  -  z"  -  R)  dR  = 


(3-36) 


This  follows  since  3p/3t  is  nonzero  at  the  boundary  of  the  patch  As.  The 
result  of  the  first  order  evaluation  of  (3-32)  is  thus 


p(r',t')  =  2  p  (z',t')  +  Y  <  p(r ' ,t' ) 


where  Y  *  t  -  2z*.  Solving  for  p  and  expanding  1/(1 -y<)  this  can  be 


.JuW.'  Wu.--  ■ .  •  . 


written 


p(r',t')  =  2  pi(2',t')  jl  +  Y<  +  y^K2  +  ...|  , 


and  making  the  substitution  (3-33) 


~  p(r'  ,t')  =  u(t  -  2z* )  |l  +  K(t  -  2z’)  +  <2(t-22')2  +  ...j  ;  (3-37) 

~  If-  U',t’>  =  6  (t  -  2z ' )  \l  +  K(t  -  2z  ’ )  +  ...j 

2  i 


+  u  ( t  -  2z ' )  <|l  +  2K(t  -  2z  ' )  +  ..  J 


(3-38) 


This  expression  for  the  surface  pressure  can  now  be  used  to  evaluate  the 
far-scattered  step  response. 


Using  (r  *n')  dS'  =  2tt  p  a  cos  0  d0  =  2n(a-zl)  d2’  ,  the  inte- 
o 

gral  (3-31)  can  be  rewritten 
t/2 

rQ  pS(t)  =  C  (a  —  z • )  6 ( _  —  , 

J0  ' 


(  —  —2 

t  -  2z‘)  ;  1  +  K (t  -  2z')  +  K 


2  ) 

(t  -  2z ’ )  +  . . . ’  dz ' 


t/2 


i 


(a  -  z' )  u(t  -  2z 


—  i  ~  —2  2  ) 

')  <1  +  2k  (t  -  2z ' )  +  3K  (t-2z‘)  +  ... 


(3-39) 


Direct  evaluation  and  the  substitution  a<  =  -1,  yields  for  the  step 
response 


r  pS  =  ai-r  -  t  (t/a)  +  4 


o  P  =  ajj  -  —  (t/a)  +  —  (t/a)  +  ...j  u(t)  . 

Differentiating  and  normalizing,  the  impulse  response  of  the  hard  sphere 
is  then 


dz’ 
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r  p  (t/a)  ,  ,  5 

h(t/a)  =  -  =  —  6  (t/a)  -  —  u(t/a)  +  —  r  (t/a)  +  . 

a  2  4  4  p 


(3-40) 


The  third  term  (ramp)  of  this  expansion  is  inaccurate  and  could  be  improved 
upon  by  the  use  of  (3-37)  and  (3-38)  in  (3-32)  as  the  next  approximation; 
and  so  on.  The  main  result  here  is  the  correction  of  the  step  coefficient 
as  given  by  the  simpler  physical  optics  theory  in  equation  (3-30) .  The 
result  of  approximating  the  impulse  response  by  the  methods  of  Section  2, 
using  a^  =  12,  is  plotted  in  Figure  3-2,  superposed  on  the  present  results. 
The  agreement,  including  the  r^  term  is  seen  to  be  excellent. 

3.2.2  Fluid  Sphere  Convex  Specular  Point 

Consider  now  the  case  where  energy  penetrates  the  target  in  the 
form  of  compression  waves.  We  then  have  two  integral  equations  (see  Sec¬ 
tion  4)  corresponding  to  pressure  outside,  p,  and  inside,  p^s 

p(r'  ,f)  =  2  p1(z'  ,t)  +  J  (R  •  n")  ^  ^  p  (r  "  ,  t "  )  ds  " 


'  Jn  f 
Jt\s 

“l 


P^r'pt1)  = 


(3-41) 


where  t"  *  t'  -  R,  t^'  =  t'  -R/c^.  The  patch  As  has  radius  y  =  t-2z*  as 
before.  The  interior  integral  is  over  the  larger  patch  As^  with  radius 
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Y  =  c^t  -  2z|.  The  above  equations  are  related  by  the  boundary  conditions 


9n  p^9n 


(3-42) 


c^,  are  che  sound  speed  and  density  inside  relative  to  outside.  The 
integrations  over  the  surface  patch  are  similar  to  that  of  the  previous 
section  and  result  in 


P(r' ,f) 


2P1 


+  YKp 


-  Y 


9n 


and 


p(r\t')  = 


IE 


C1  P  +  YP1C1  9n 


where  the  last  equation  is  obtained  from  the  last  of  (3-41)  and  the  appli¬ 
cation  of  the  boundary  conditions.  Again,  we  use  a  step  function  as  the 
incident  field.  Solving,  we  obtain  on  the  surface 


P  C 

(r\f)  =  2^  +1p-1-c-  j  (l  +  e  +  e2  +  ...)  u (t  -  2z ' ) 


*2(r r^)  (> 


+  e  +  e  + 


(t  -  2z ' )  , 


(3-43) 


where  we  have  defined 


-  (  pi  —  1 


<  c1  Y 


The  integral  equation  for  far  field  scattering  from  an  elastic  or  fluid 


target  is 


r  pS(t) 
o 


8p(r' ,t‘) 

at' 


dS' 


4TT 


f 


8p(r',f) 


3n 1 


dS' 


(3-4< 


Substituting  (3-43)  into  this,  results  in  terms  similar  to  that  for  the 
hard  target  step  response  except  for  constants.  The  result  is  that  the 
impulse  response  of  a  smooth  curved  specular  point  for  a  fluid  target  is 
given  by 


h(t/a) 


(3-4E 


where 


6  = 


pici 


Vi + 1 


Note  that  as  P1  “  (3-45)  reduces  to  the  sound-hard  case.  Also,  the 
factor  (P^c^  -  ll/(p^c^  + 1) ,  by  which  the  impulse  term  differs  from  the 
sound-hard  (perfectly  reflecting)  case,  is  the  same  as  the  reflection 
coefficient  found  in  (3-27)  for  normal  reflection  from  an  elastic  plane 
surface. 


3.2.3  Extensions  of  These  Techniques 

At  the  next  level  of  complexity  it  is  possible  to  determine  by 
these  analytic  methods  the  response  of  an  elastic  specular  point.  The 
expressions  become  very  complicated  however  when  shear  wave  effects  are 
included;  hence  this  approach  was  not  pursued  here.  A  direct  approach  is 


the  numerical  evaluation  of  the  integrals  of  the  preceding  section  usj.ng 
a  very  short  smoothed  impulse  incident  pressure.  The  result  is  the  leading 
edge  of  the  impulse  response  and  requires  the  numerical  integration  over 
only  a  small  part  of  the  target  surrounding  the  specular  point.  Further¬ 
more,  the  technique  may  be  applied  to  a  bistatic  configuration  and  also 
to  refraction.  This  approach,  known  as  the  local  influence  theory,  is 
based  on  the  observation  that  at  a  given  time  only  a  few  identifiable 
points  are  active  as  scatterers;  see  Figure  3-3.  (The  approach  has  been 
used  successfully  for  the  electromagnetic  case  [10].)  Attempts  to  imple¬ 
ment  this  technique  for  elastic  targets  have  been,  to  date,  thwarted  by 
the  numerical  difficulties  further  described  in  Section  4.  It  is  felt 
that,  once  these  difficulties  are  removed,  these  calculations  of  scattering 
center  responses,  combined  with  the  glory-wave  model  to  be  described,  will 
yield  a  fully  quantitative  but  greatly  simplified  scattering  model. 

3.3  GLORY -WAVE  EFFECTS 

The  glory  phenomenon  is  the  halo  surrounding  the  viewer's  own 
shadow  cast  by  the  sun  upon  a  mist  (or  a  cloud,  when  observed  by  airplane) . 
It  is  caused  by  light  entering  water  droplets  and,  after  multiple  internal 
reflections,  being  returned  to  the  observer  [11][12].  It  is  a  monostatic 
effect  as  distinguished  from  the  rainbow  effect  which,  although  also 
involving  internal  reflections,  is  bistatic,  resulting  from  local  maxima 
in  scattering  intensity  at  certain  bistatic  angles.  The  phenomena  involve 
surface  (creep)  waves  as  well  as  Snell's  law  reflection  and  refraction. 

The  nomenclature  is  here  applied  to  the  acoustic  case.  (Since 
the  refractive  index  for  most  acoustical  situations  will  be  greater  than 
1  while  for  the  optical  case  it  is  generally  less  than  1,  there  are  some 
qualitative  differences  between  the  two  cases.)  A  semiquantitative  under¬ 
standing  of  the  phenomenon  can  be  obtained  using  ray-optics  theory. 

Consider  the  glory  path  as  illustrated  in  Figure  3-4  for  a 
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FIG.  3-3  Local  influence  theory:  reflection  and  refraction, 
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sphere.  At  each  change  in  direction  Snell's  law  is  obeyed 


sin  8,  sin  8„  sin  0, 
12  3 


(3-46) 


where  c^,  c^,  c^  are  the  compression,  shear,  and  external  velocities 
respectively.  Each  chord  of  the  path  is  thus  the  base  of  an  isosceles 
triangle.  The  requirement  is  that  after  k  traverses  around  the  sphere 
the  rays  return  to  the  source.  Let  and  m2  be  the  number  of  chords  at 
c^  and  c2  respectively.  Then 

203  +  m1('rT  "  201)  +  n>2(iT  -  282j  =  2TTk  (3-47) 


With  the  constraints  (3-46) ,  this  equation  has  at  most  one  solution  8^ 
for  each  set  (m^,  m2 ,  k) . 

The  extra  time  required  to  traverse  this  path  compared  to  direct 
reflection  from  the  center  of  the  sphere  is 


t  ,  =  2a  I  —  cos  0,  +  —  cos  0  -  cos  0 

n'lm2k  \  C1  1  c2  2  3 


(3-48) 


where  a  is  the  sphere  radius. 


The  cases  m^  or  m2  -*■  00  are  identified  as  creep  waves  along  the 

(inner)  surface  of  the  sphere,  which  are  launched  when  8^  is  at  one  of 

the  critical  angles.  By  this  model,  then,  the  impulse  response  of  an 

elastic  sphere  is  a  sequence  of  pulses  occurring  at  the  times  t  .  . 

#  ni^  f k 

The  magnitudes  and  shapes  of  these  pulses  depend  upon  the  detailed 
interaction  of  the  local  wavefront  with  the  boundary  at  each  of  the  re¬ 
flection  points. 


3.4  GLORY  MODEL  EVALUATION 

Consider  first  the  fluid  sphere  response  (that  is,  considering 
only  compression  waves).  Figure  3.5  illustrates  two  cases,  comparing  the 
smoothed  impulse  response  with  the  impulses  predicted  by  the  glory  wave 
model.  It  is  seen  that  the  agreement  is  virtually  perfect.  All  possible 
glory  paths  are  observed  with  decreasing  amplitude  in  proportion  to  the 
number  of  reflections. 

In  Figures  3.6  and  3.7  the  model  is  compared  with  the  elastic 
target  response.  First,  note  that  the  pure  compression  modes  as  observed 
for  the  fluid  target  are  still  observed  (for  the  sake  of  clarity  not  all 
of  them  are  indicated) ,  though  these  effects  are  overwhelmed  by  the  shear 
effects.  The  agreement  between  the  glory  modes  and  the  main  features  of 
the  response  is  only  approximate.  It  is  significant  that  the  change  in 
times  of  occurrence  with  changing  c 2  and  c^  agree. 

It  is  seen  from  the  figures  that  the  dominant  modes  are  those 
with  m^  *  1,  m2  =  2'  or  m2  =  °°  (creeP) •  There  is  little  time  difference 
between  (ir^,  m^)  =  (2,1)  and  (°°,1)  or  between  (1,2)  and  (°°,2)  ,  where 
m^  =  00  means  that  the  path  at  speed  c^  is  an  internal  creep  path.  A  more 
quantitative  theory  is  required  to  determine  the  relative  contributions 
from  those  modes.  These  principal  modes  are  illustrated  in  Figure  3-8. 

3.5  SIMPLE  MODEL 

It  is  apparent  that  pure  compressional  modes  can  be  neglected 
in  a  single  scattering  model  for  elastic  solids.  We  will  base  our  model 
on  the  paths  a,  b,  c,  d  of  Figure  3-8.  For  path  (a) ,  the  specular 
response,  the  return  occurs  at  time 


t 

a 


(3-49) 
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FIG.  3-5  Glory  wave  model  —  fluid  sphere  response. 


(a)  SPECULAR 


FIG.  3-8  Principal  scattering  modes  involving  shear  waves. 


For  paths  (b)  and  (c)  we  are  at  the  critical  angle  for  c  .  That  is 

0,  =  90°  and  8.  ,  =  sin  1  (c  /c,)  and  0  =  sin  1  (c,/c.) .  The  creep  path 

1  3  c  1  31  2  21 

length  travelled  at  c^  is  just  sufficient  to  satisfy  the  requirement  that 

the  ray  return  to  the  source.  (Actually,  the  creep  wave  at  c^  is  continually 

radiating  so  that  all  orientations  of  the  path  at  c^  occur.  Also,  at  other 

observation  angles,  the  same  model  can  be  expected  to  apply  with  a  different 

arc  length  at  c^ . )  The  arc  at  spans  in  case  (b)  :  =  8^  +  (tt/2  -  0^  ) 

and  in  case  ,c) :  0,  =  20_  -  0_  ,  . 

1  2  3cl 

The  return  times  are  thus 


t  _  2a  !l  +  C°S  92  _  C°S  63cl 
b  3  C1  C2  C3 


t  =  2a  +  2  C°S  62  _  C°S  93cl 

C  C1  C2  C3 


(3-50) 


(3-51) 


Path  (d)  has  the  critical  angle  associated  with  c2»  ®3c2  =  s^-n  ( ^  ' 

and  half  arc  length  Y„  =  ( tt  —  0,  .) .  Hence  the  return  time  is 

2  3c2 


t,  =  2a  -1 
d  c 


X  .  COS  63c2 


(3-52) 


A  computer  program  was  written  to  implement  this  model.  The 
response  for  each  mode  is  assumed  to  be  of  the  form  illustrated  in 
Figure  3-9,  namely  impulse,  step,  ramp,  and  backswing  (P,  S,  R,  and  E) 
such  that  the  total  area  (DC  response)  is  zero. 


:,tl  Pk  6(t-tk)  *  Sk  u(t-tk)  .  \(t-tk)[u(t-tk)  -  u(ftk.At)j 


*  Ek  u(t_tk  *  At) 


(3-53) 
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At  this  stage  P,  S,  and  E  are  inputs  to  the  program;  with  At  =  -2P/(E+S), 
R  =  (E  -  S) /At .  It  is  expected  that  development  of  the  methods  in  Sections 
3-1  and  3-2  will  yield  theoretical  values.  The  times  of  occurrence  of 
these  pulses  are  computed  from  the  relations  (3-49)  to  (3-52) .  It  is  seen 
in  Figure  3-10  for  aluminum  that  a  simple  representation  for  the  impulse 
response  can  be  found  in  this  way.  In  the  frequency  domain  (Figure  3-11) , 
the  model  and  actual  response  (as  computed  by  Section  2  methods)  agree 
in  the  main  features.  Brass  is  illustrated  in  Figures  3-12  and  3-13. 


SECTION  4 

SPACE-TIME  INTEGRAL  EQUATION  APPROACH  (SOLID  TARGETS) 


A  major  advantage  of  an  integral  equation  formulation  of  the 
scattering  problem  is  that  it  yields  exact  responses  for  targets  of  arbi¬ 
trary  three  dimensional  shape.  The  integral  equation  technique  (STIE) 
used  at  the  Sperry  Research  Center  operates  in  the  time  domain,  solving 
the  problem  by  stepping  in  time  rather  than  by  matrix  inversion.  Using 
the  excitation  described  in  Section  2.2,  the  result  is  a  smoothed  impulse 
response.  Such  a  time  domain  response  is  highly  suggestive  of  target 
shape.  If  desired,  the  response  can  be  deconvolved  with  the  incident 
pulse  and  transformed  to  the  frequency  domain  to  obtain  a  frequency 
response.  The  latter  is  valid  for  values  of  ka  from  zero  to  an  upper 
limit  determined  by  the  frequency  content  of  the  excitation  pulse.  The 
latter  is  limited  primarily  by  computation  time.  To  date  most  experience 
with  this  type  of  solution  has  been  with  electromagnetic  scattering  and 
with  acoustic  scattering  from  sound  hard  and  sound  soft  targets  [6,  13,  14]. 
It  has  been  found  that  practical  values  on  the  upper  frequency  limit  are 
ka  =  6  to  12. 

To  date  the  STIE  approach  has  been  applied  with  complete  success 
to  impenetrable  targets  such  as  conducting  solids  and  thin  surfaces  in  the 
electromagnetic  case  and  to  sound  hard  and  sound  soft  targets  in  the 
acoustic  case.  Results  have  been  published  [6,  13,  14]  for  three  dimen¬ 
sional  targets  of  varying  complexity. 

In  the  present  effort,  this  technique  is  extended  to  penetrable 
targets.  The  approach  is  reported  in  detail  in  reference  1,  and  will  here 
be  outlined  in  a  somewhat  more  convenient  form.  At  time  of  writing,  tar¬ 
gets  penetrated  only  by  compression  waves  are  treated  successfully,  whereas 
numerical  difficulties  plague  the  solution  when  shear  waves  are  added. 

We  will  first  concentrate  on  the  fluid  target  solution  and  present 
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computational  results. 


4.1  FLUID  TARGET  FORMULATION 


In  an  ideal  fluid  the  excess  pressure  satisfies  the  wave  equa¬ 


tion 


2  „2 

P  =  c ,  v  p  , 


(4-1) 


where  c  is  the  speed  of  sound.  The  pressure  is  related  to  particle  dis- 
■*  -*• 

placement  u  by 


Vp  =  -p3  u 


(4-2) 


where  is  the  density  (see  [1],  Section  2.1) 


4.1.1  Derivation  of  Integral  Equation 


Consider  a  closed  surface  S  in  the  medium  with  pressure  sources 
p1,  originating  outside  of  S.  Using  the  Green's  function 

G(r|r\t|t)  -  j  -  R/c3  -  tj 

where  R  =  r  -  r' ,  a  solution  to  the  wave  equation  can  be  developed,  known 
as  the  Kirchhoff  integral  representation: 

PU,t>  =  pXU,t)  -  ~  f  -9£-fn-r-T-  -  ("•  •  R)  l  pU'  , T) |  ds 1  (4-3) 


where  the  integral  is  over  the  surface 

n'  is  the  unit  normal  out  of  the  surface 
R  is  the  magnitude  of  R  ■  r  -  r' 
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R  is  the  unit  direction  of  R 


r  _  1  1  j 

-  r2  Rc3  9l 


T  =  t  -  R/c, 


The  above  is  valid  for  r  outside  of  S. 


Let  the  inside  of  S  be  source-free  and  filled  with  a  fluid  with 
properties  c  ,  p  ;  then  (maintaining  the  same  direction  for  n) 

-7?  jf  { 51^  •  *>  h 


for  r  inside  S ,  where 

l  -L,  J— i 

1  B2  SC1  9T1 
Tj_  ■  t  -  R/C 

and  p^,  9p^/9n  are  the  values  just  on  the  inside  of  the  surface. 

That  is,  the  pressure  can  be  calculated  everywhere  if  the  values 
of  p  and  9p/9n  are  known  on  S.  In  particular,  in  the  far  field  as 
p  -+  r  -*•  oo,  the  scattered  part  of  the  field  is  given  by 


r  p  (r, t)  = 


■i  /{- 


dp  ,~+ .  /v  /v  dp  -*• 

aj  (r  »T)  +  (n'  *r)  (r 


'  ,T)jdS’ 


To  evaluate  this  expression,  it  is  first  required  to  solve  (4-3)  and  (4-4) 
on  S.  A  limiting  procedure  must  be  applied  to  make  these  valid  on  the 
surface  since  the  term  (n  •  R)  L p in  the  integral  becomes  singular  as  r’  +  r 
This  singularity  is  removed  in  the  following  way:  consider  the  singular 
term  separately  and  define  the  integral  written  as  J  dS'  to  have  its 


Now  require 


singular  term  removed.  (Usually  this  is  written  as  f  dS'.) 

■+  *  -> 

that  (4-3) ,  the  "outside  equation"  results  in  p(r)  =  0  for  r  inside  S. 

This  requirement  remains  to  be  justified.  Comparing  the  results  of  eval- 
uating  (4-3)  for  r  just  inside  and  just  outside  of  S,  we  see  that  the  non¬ 
singular  parts  of  the  equation  I  _  give  the  same  value  at  these  two  points, 

NS 

but  that  the  singular  portion  changes  sign  (due  to  r  •  n) .  That  is,  we 
have 


r  outside:  p  =  p1  +  I  -  Ig 
x  inside:  0  =  p  =  p1  +  INS  +  ls 


(4-6) 


Combining  these,  we  obtain  p  =  2P1  +  2  I  .  The  singular  point  has  been 

-*■ 

removed  and  we  have  introduced  a  factor  of  2 .  We  thus  obtain  for  r  on  the 
(outside  layer  of  the)  surface  S: 


p(r,t)  =  2p  (r,t) 


_  J_  f  f±Jz. 

2tt  Js  \  R  9n'  lr 


-  (n'  •  R)  L  p(r,T) 


|  dS' 


(4-7) 


This  equation  contains  two  unknowns  p  and  9p/9n  so  that  another  relation 
is  needed  to  solve.  In  the  case  of  an  impenetrable  target  the  requirement 
p  =  0  inside  is  exactly  physically  true  so  that  p  and  9p/9n  in  (4-7)  repre¬ 
sent  actual  physical  quantities.  The  additional  boundary  condition 
9p/9n  =  0  for  sound-hard  or  p  =  0  for  sound-soft  then  permits  solution  of 
(4-7)  . 


In  the  case  of  the  fluid  target  we  apply  an  analogous  argument, 

*4  4 

requiring  inside  equation  (4-5)  to  yield  p(r)  =  0  for  r  outside,  resulting 
in 


s  '  '  (4-8) 
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where  K.  is  the  average  curvature  of  the  patch  at  r..  The  derivation  of 
—x  x 

(4-11)  is  straightforward,  by  approximating  the  patch  by  a  circle  with 
radius  *y  .  For  (4-12)  it  is  necessary  to  expand  R  and  n  in  curvilinear 
coordinates  (see  Appendix  in  reference  13) .  A  similar  numerical  expres¬ 
sion  is  written  for  (4-8) .  Writing  1  and  I  for  these  sums  (the  "non-self 
integrals")  and  applying  the  boundary  conditions,  we  obtain  the  system 


(i  -  £.)  p(vtj)  +  Yi  t  (vh) '  2pl(  vh)  ■ !( Vh) 


(4-13) 


This  can  solved  simultaneously  for  p,  3p/3n  for  time  t_.  at  each  r^  on  S. 
Note  that  the  non-self  integrals  can  be  evaluated  at  time  t_.  since  the 
arguments  T  are  earlier  than  t_.  .  Hence  (4-13)  can  be  solved  for  all  time 
by  stepping  in  time. 


The  stepping  in  time  procedure  is  valid  so  long  as  the  distance 
Ar  to  the  nearest  patch  is  greater  than  c^At  (or  c^At,  whichever  is  greater) 
for  At  =  time  step.  It  may  be  desirable  to  use  smaller  patches  in  regions 
of  rapid  change,  for  which  Ar  is  smaller  than  the  above  limit.  For  these 
cases  a  local  iterative  procedure  is  used  involving  only  these  nearby 
patches:  first  the  contributions  from  all  the  non-near  (properly  retarded) 
patches  are  calculated  to  each  patch.  These  values  are  then  used  to  deter¬ 
mine  the  near-patch  contributions  by  simple  linear  interpolation  between 
the  new  values  and  those  at  the  previous  time  point.  Then,  one  more 
iteration  of  this  process  is  sufficient  to  assure  convergence.  This  pro¬ 
cedure,  involving  only  the  nearby  patches  is  quick  and  does  not  add  much 
time  to  the  overall  solution. 


Typical  values  for  a  target  dimension  is  a  =  1  (representing, 

for  instance,  cylinder  radius)  and  for  the  pulse  length  parameter  a  =  1 

n 

(see  equation  2-15).  Then,  for  time  expressed  in  sound-meters  (c^  =  1) , 
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a  time  step  of  At  =  .2  (or  even  .3)  is  adequate  to  obtain  good  results  in 
the  case  of  sound-hard  or  soft  targets.  Since  the  pulsewidth  measured 
between  2%  points  is  about  W  =  4/a^  =  4,  this  represents  about  20  time 
steps  per  pulsewidth.  The  typical  patch  separation  then  is  somewhat 
larger  than  0.2.  For  the  fluid  case,  the  time  step  must  be  reduced  by 
about  the  ratio  c.j/c^.  It  was  found  that  At  =  .075  was  sufficiently  fine. 

In  the  integral  (4-10)  it  is  necessary  to  interpolate  in  the 
table  of  previously  computed  values  of  p  to  obtain  p  and  3p/3t  at  time  T. 
One  way  to  do  this  is  by  fitting  a  fourth-order  polynomial  to  the  five 
values  surrounding  t.  This  worked  very  successfully  and  accurately  for 
the  sound-hard  case.  However,  numerical  inaccuracies  tend  to  develop 
toward  the  end  of  the  response  in  time  which  may  be  amplified  by  high 
order  polynomial  fitting,  resulting  in  a  rapidly  growing  instability. 

This  effect  was  particularly  noticeable  for  the  sound-soft  case.  This 
problem  can  be  eliminated  by  linearly  smoothing  the  five  points  and  then 
fitting  a  simple  quadratic  to  the  inner  three  points.  This  method  causes 
a  small  loss  in  accuracy  and  hence  requires  the  use  of  a  somewhat  smaller 
At  (say  .1  for  sound  soft,  or  .05  for  fluid  with  c^/c^  =  4) . 

4.1.3  Responses  of  Right  Circular  Cylinders 

Responses  were  computed  by  these  methods  for  a  fluid  right  cir¬ 
cular  cylinder  with  the  properties  of  aluminum  and  brass  but  neglecting 
the  shear  velocity: 


P1 

c 

"Fluid 

Al" 

2.70  g/cm3 

6200 

"Fluid 

Brass" 

8.39 

4400 

Medium 

(Water) 

1.00 

1480 

the  incident  pulse  was  a  *  1.  Results  are  plotted  in  Figures  4-1  through 
4-3.  The  cylinder  has  length/radius  ratio  L/ a  =  6/1.  Time  is  in  units  of 
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c^t/a.  Consider  the  0  aspect;  the  initial  part  of  the  response  is  from 
the  flat  end  of  the  cylinder  and  looks  somewhat  like  the  derivative  of 
the  incident  pulse,  centered  on  t  =  -6.  The  small  peak  at  about  t  =  -4 
can  be  interpreted  as  interation  with  the  edge.  This  part  of  the  response 
is  similar  in  shape  to  the  response  of  a  hard  cylinder  (see  [13]  or  [1]) 
except  for  a  reduction  in  amplitude.  The  small  peak  at  t  =  -4  is  also 
somewhat  earlier,  indicating  that  edge  effects  occur  at  both  c^  and  c^. 

The  next  major  part  of  the  return  is  the  creep  return  at  t  =  6  to  8.  This 
has  been  changed  in  shape  and  moved  slightly  earlier  ir.  comparison  with 
the  hard  cylinder  response.  In  between  these  two  extreme  times,  the  hard 
cylinder  has  a  near-zero  response,  whereas  the  fluid  target  exhibits 
ringing.  The  latter  must  be  interpreted  as  the  result  of  internal  reflec¬ 
tions  (end  to  end)  of  compression  waves  at  c^.  As  the  aspect  angle  changes 
the  two  parts  of  the  response  move  closer  together. 

A  similar  set  of  responses  is  given  for  brass  in  Figure  4-2. 

Also  given,  in  Figure  4-3,  is  a  comparison  for  each  of  these  targets  at 
the  0°  and  4°  aspects.  The  difference  is  extremely  small. 

4.2  ELASTIC  TARGET  FORMULATION 

A  complete  elastic  target  formulation  must  also  include  the 
effect  of  internal  shear  waves.  As  seen  in  Section  2,  these  can  be  the 
dominant  effects.  This  formulation  is  treated  in  detail  in  [1]  and  will 
only  be  cursorily  treated  here.  Instead  of  pressure  it  is  easier  to  deal 
with  the  divergence  of  displacement.  The  two  quantities  are  proportional: 


P1  =  ”  C1  Pl(7  ’  u) 


The  equation  of  motion  inside  the  solid  is 


u  =  c^  V(V  •  u)  -  c^  V  x  V  x  u  , 


(4-14) 


(4-15) 


vh".:nA" 


from  which  it  may  be  found  (by  taking  the  V  x  and  V  •  of  the  above)  that 
both  (V  •  u)  and  (Vxu)  separately  satisfy  wave  equations: 

•  • 

(v  x  u  )  =  c2  V2(Vxu)  (4-16) 

-f 

In  addition  to  equations  (4-7)  and  (4-8) ,  with  replaced  by  (V  •  u) ,  we 
have  a  vector  Kirchhoff  integral: 


Vxu  (r,t)  = 


,J-  f  fi_L 

2 7T  J  JR  3n‘ 


(Vxu’)  -  (n'  •  R)  L  (V  x 


u’)| 


dS'  (4-17) 


for  r  on  S.  Here,  the  equation  has  already  been  specialized  to  the  boun¬ 
dary  S. 


Vxu'  means  V  x 


T 


2 


t  -  R/c2  . 


Another  form  for  this  equation,  which  may  be  useful  in  practice  is 
Vx”  =  -JH  f  h  (n'xVxVxu')  +  L2  |(fi>  .Vxu- 

Q  ' 


)R+  (n'xVxu’)xRjjd 
(4-18) 


This  form  is  obtained  from  (4-17)  by  vector  manipulations  and  the  use  of 

the  divergence  theorem.  (See  Jackson  [8],  pp.  283.)  It  may  be  useful  to 

observe  that  another  equation  can  be  written,  namely  for  the  (VxVxu). 

The  quantities  (Vxu),  and  (VxVxu)  now  play  somewhat  similar  roles  to 
-*•  -*■ 

E  and  H  of  electromagnetics,  so  that  techniques  found  useful  there  are 
possibly  applicable  here. 


The  boundary  conditions  are  the  two  stated  in  Section  4.1,  con¬ 
tinuity  of  normal  stress  and  normal  displacement,  although  the  first  looks 


somewhat  more  complicated  for  elastic  solids  (see  [1]) .  An  additional 
boundary  condition  on  the  shear  states  that  the  tangential  shear  stress 
is  continuous  across  S.  Since  this  quantity  is  zero  in  the  fluid,  it  is 
zero  at  the  solid  side  of  S.  The  boundary  conditions  may  be  written  (see 
[1]  for  derivation;  here  the  form  in  which  these  relations  are  written  is 
more  precise) : 


-P  =  P^i  ~  2°2 1  * 


0  3u 
2  n 

+  2  °1  C2  IS" 


<4-19) 


3p  2  3  2  •  • 

— 2  =  c  —  (V  •  u)  -  c_  (V  x  V  x  u)  =  u 

P33n  1  3n  2  n  n 


(4-20) 


0  =  (V  x  u)  -  2  —  (nxu) 
tan  dn 


(4-21) 


where  n  is  the  outward  normal  direction.  The  set  of  2  scalar  and  one  vec¬ 
tor  STIE's  plus  the  two  scalar  and  one  vector  boundary  equations  can  be 
solved  for  the  unknowns  p,  V  •  u,  Vxu  and  their  normal  derivatives.  How¬ 
ever  it  requires  the  formation,  numerically,  of  certain  space-derivatives 
on  the  surface  of  S ,  and  the  numerical  solution  of  the  equation  of  motion 
(4-15)  [1].  The  need  for  differentiation  on  a  curved  surface  results  in 
numerical  instabilities  in  the  solution,  which  have  to  date  not  been  solved. 
Only  the  initial  portion  of  the  response  from  the  region  around  the  specu¬ 
lar  point  can  be  computed  successfully. 


Figure  4-4  illustrates  the  status  of  the  space-time  integral 
equation  calculation.  Plotted  is  the  surface  pressure  as  a  function  of 
time  pfc^t/a)  for  the  smooth  impulse  incidence  at  several  points  on  the 
surface.  The  pressure  as  predicted  by  the  classical  solution  of  Section  2 
is  plotted  in  part  (a) ;  the  current  status  of  the  integral  equation  solu¬ 
tion  in  part  (b) .  The  shape  of  p  more  or  less  follows  that  of  the  incident 
pulse  with  significant  deviations.  In  particular,  at  the  backside 
(8  =  180°)  there  is  a  bump  a  little  before  t  =  0  which  represents 
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transmission  through  the  target  at  the  shear  velocity,  c^.  It  is  seen  that 
the  initial  part  of  the  STIE  calculation  is  correct;  in  particular  the 
bump  at  0  =  180°  is  computed  correctly-  The  latter  is  the  respect  in  which 
elastic  targets  are  distinguished  from  hard  or  fluid  targets.  Unfortunately, 
inaccuracies  which  rapidly  develop  into  an  instability  start  near  0  =  45°. 
These  inaccuracies  have  been  identified  as  numerical  differentiation  prob- 
lems;  particularly  since  the  quantity  differentiated,  u,  becomes  large  as 
time  goes  on,  although  the  derivatives  to  be  found  remain  relatively  small. 
This  problem  has  not  yet  been  circumvented.  The  solution  must  lie  in  a 
formulation  which  avoids  space-differentiation  of  the  displacement  u.  The 
far  field  computed  from  these  surface  values  is  shown  in  Figure  4-5,  with 
the  classical  result  shown  dashed. 


4.3  THIN  SHELL  SPACE-TIME  INTEGRAL  EQUATION  FORMULATION 


A  direct  approach  to  solving  the  problem  of  scattering  from  an 
elastic  shell  with  fluid  interior  and  exterior  is  to  formulate  integral 
equations  over  both  boundaries  and  applying  boundary  conditions  on  each. 
This  two-boundary  approach  is  taken  in  the  classical  solution  of  the  next 
section.  A  computationally  simpler  approach  is  to  attempt  to  find  the 
limiting  boundary  condition  between  inside  and  outside  media,  and  avoiding 
the  need  for  solving  any  equations  in  the  shell  itself. 


Consider  the  sphere  of  Figure  4-6  with  outer  and  inner  radii 
a  and  b.  The  boundary  conditions  can  be  written  as 

a 


-p  -  C^c2  -  2c2)  (V-I)3  *  2  Pj  <=1  (-5^) 
-P„  '  P,(c!  '  2c2  )  IV‘“>b  *  2  P1  c2  (-37) 


(4-22) 


and 
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••a 

u 

r 


Sp  2  j_ 
p^9r  ~  C1  9r 


(V  •  u)a 


c?  (V  x  V  x  u) a  , 
2  r 


•  *b 

u 

r 


3p4  _  2  _9_ 

p^3r  ~  C1  9r 


(V  •  u)b 


cl  (V  x  V  x  u)b  ; 
2  r 


(4-23) 


where  the  superscripts  a  and  b  refer  to  the  outer  and  inner  surfaces.  If 
we  write  (VxVxu)^  and  the  boundary  condition  (4-21)  in  spherical  coor¬ 
dinates  ,  we  obtain 


(VxVxu)a  =  2  ~  (V*u)a  -  2  ~  |  -77 


r  [7^  (r\) 


(4-24) 


We  use  this  to  rewrite  (4-23)  : 


••a 

u 


P39r 


-  (4  -  -22) 


[7  £  (r2“t)] 


3p, 


P43 


1  /  2  _  2\  9  -*■  b  29  fl  9/2  \"| 

x  "  \C1  ‘  2C2)  97  (V  ‘  U)  +2C2^  (r  Ur)J 


(4-25) 


9P„ 


3  "  r4 

The  aim  is  to  find  relations  between  p,  7^7  p4,  -7^-  by  applying  a  limiting 
process  so  that  the  problem  reduces  to  essentially  the  fluid-in-fluid 
problem  of  Section  4.1,  but  with  different  boundary  conditions. 

First,  suppose  the  shell  thickness,  h,  sufficiently  thin  so  that 
we  can  use  a  first  order  Taylor  expansion;  for  example 


(V  •  u)  3  =  (V*u)b  +  h  (V*u) 


(4-26) 


We  will  then  rewrite  the  boundary  conditions 
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P4  -  P  -  <V>(cJ  -  2c2)  £  (V-:>  *  2  Plh 


a2 

0  a  u 
2  _ r 

'2  9r2 


3p4 

p43r 


_§E_  = 

P39r 


'K  '  2C2>^2  *  2hC2  ^2  [7^  (“r)] 


(4-27) 


(4-28) 


If  both  h  and  p  are  reduced  to  zero  the  problem  reduces  to  the  fluid-in¬ 
fluid  problem  with  no  shell.  As  a  next  level  of  approximation,  let 
b  a 

h  -*•  0  and  let  u  =  u  .  Then 
r  r 


P4  ‘ 


’  0lhK  -  24)ir  IV 


u) 


u 

r 


3p4  5p  _  /  2 

p43r  '  p33r  =  \C1 


(V  •  u)  . 


(4-29) 


(4-30) 


or 


p4  -  p  =  pl  h  ur  •  (4-29') 

This  represents  the  equation  of  motion  of  a  membrane  with  no  internal 
restoring  forces.  We  can  add  restoring  forces  due  to  an  ambient  pressure 
difference,  Ap,  between  inside  and  outside  which  results  in  a  tension 
T  *  Ap  (a/2) .  Then  we  have 

•  •  2 

p  -p=p  hu  -  T  V  u  .  (4-29") 

4  1  r  s  r 

This  is  the  equation  of  motion  of  a  membrane  under  tension  and  can  be 
solved  in  conjunction  with  the  fluid-in-fluid  integral  equations. 

We  have  simplified  more  than  desired,  however,  since  we  wish  the 
elastic  properties  of  the  shell  to  be  included.  Returning  to  equations 
(4-23) ,  we  note  that  we  can  also  write  by  virtue  of  the  boundary  conditions 
Cl] 
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(4-31) 


(V  x  V  x  u)  =  -  2  V2  u  , 
r  s  r 


2 

but  it  must  be  realized  that  the  meaning  of  is  dependent  on  the  coor¬ 

dinate  system.  In  spherical  coordinates  (with  axial  symmetry) : 


s 


u 

r 


.2 
3  u 


3u, 


2  2  2 
r  38  r  38 


cos  8 

2  ■  « 
Rin  n 


(K 

\  36 


(4-32) 


the  above  is  actually  a  definition  of  V  ,  obtained  from  writing  out 
(VxVxu)  and  applying  the  boundary  condition  in  spherical  coordinates. 
We  then  write  (4-23)  as 


•  •  2  3  ■+■  o  ? 

u  =  c,  -r—  (V  •  u)  +  2  c_  V  u 

r  1  3r  2  s  r 


(4-33) 


and  use  this  to  eliminate  3/3r  (V  *  u)  in  equation  (4-29) .  The  result  is 


P,  -  P  =  (P,h)  (1  -  2  K  -  2  -|  V2sur 


(4-34) 


•  •  3p 

ii  =  -  —  - —  = 


9P, 


P33r  p43r 


(4-25') 


These  equations  can  be  used  as  the  boundary  conditions  with  the  system  of 

two  space-time  integral  equations  as  written  for  the  fluid  case  of  Section 

•  • 

4.1.  It  requires  the  time-integration  of  u^  (which  is  no  problem)  and  the 
numerical  differentiation  in  the  surface  of  u^.  Unfortunately,  as  was  the 
case  for  the  elastic  target,  the  numerical  differentiation  causes  insta¬ 
bilities. 


In  Figure  4.7  is  plotted  the  result  of  the  numerical  implementa¬ 
tion  of  the  space-time  integral  equation  solution  using  the  membrane 
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boundary  condition  between  similar  fluids.  As  will  be  seen  in  the  next 
section,  this  result  is  not  a  good  model  for  a  thin  shell.  While  the 
general  nature  of  the  response  is  similar,  its  dominant  frequency  differs 
by  a  factor  of  2.  It  is  therefore  necessary  to  solve  the  problem  remaining 
with  implementing  boundary  condition  (4-34) . 
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SECTION  5 


SPHERICAL  HULLED  TARGETS 


5.1  CLASSICAL  SOLUTION  FOR  THICK  SHELL 

The  fluid-filled  spherical  thick  shell  problem  can  be  solved  by 

expansion  in  terms  of  eigenfunctions  in  a  way  completely  analogous  to  the 

elastic  solid  of  Section  2.1  [15].  Let  the  outer  and  inner  radii  be  a,b; 

and  let  the  core  have  constants  c . ,  p  .  In  addition  to  the  incident  and 

4  4 

scattered  pressure  expansions,  we  have  a  compressional  wave  in  the  sphere 

core  (p  ,  c .) .  The  eigenfunctions  must  be  bounded  at  the  origin  and  are 
4  4 

therefore  j  (kr)  P  (cos  0).  In  the  shell,  on  the  other  hand,  there  no 
n  n 

longer  is  the  requirement  that  the  eigenfunction  be  bounded  at  the  origin. 
Or,  alternatively  we  can  say  that  the  shell  exhibits  both  inward  and  out¬ 
ward  traveling  waves.  The  shell  solution  thus  contains  both  j  and  n 
.  ’  n  n 

(or,  both  n  and  n  ,  whichever  is  more  convenient) . 
n  n 

-iiot 

Suppressing  the  time  dependence  e  ,  we  thus  have 


. 

p1  (r)  =  P  (2n  +  l)(-i)n  j  (kr)  P  (cos  0)  (5-1) 

o  n  n 


n=0 


t  i  a 


pS(r)  =  P  \ '  c  h^  (kr)  P  (cos  0) 
o  n  n  n 


(5-2) 


n=0 


r  5  a 


p4(^  =  PoEfn  jn(V)  Pn<COS  6) 

n=0  r  <  b 


(5-3) 


=  S  ("n  jnM  +  dn  nn(V))  V 


n=0 


cos  0)  (5-4) 

b  <  r  <  a 
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tiMiffiHif  nrr  iV  fi  n  i  ■ 


A.  (r) 


L 

n=0 


(b»  Jn(k2r)  +  en  "„(V)) 


dP 
_ n 

d6 


(5-5) 


b  <  r  s  a 


Again  we  have  that  in  the  elastic  shell 

u(r)  =  -  Vtj;  +  V  x  A  b  <  r  <  a 


(5-6) 


p  u  (r)  =  (X  +  2y)  V(V  •  u)  -  y  VxVxu  b  <  r  £  a  (5-7) 

and  that  the  total  pressure  in  the  external  medium  is 

p  =  p1  +  pS  r  £  a  (5-8) 

The  unknowns  are  a  ,  b  ,  c  ,  d  ,  e  ,  f  for  each  mode  n.  The  boundary 
n  n  n  n  n  n 

conditions  (2-6) ,  (2-7)  and  (2-8)  are  applied  both  at  r  =  a  and  at  r  =  b. 
The  result  is  a  system  of  6  equations  in  6  unknowns  analogous  to  the 
system  (2-10) : 


Mnzn  =  Vn 


(5-9) 


where  Z  is  written  for  the  vector  of  unknowns  (a  ,  b  ,  . f).  The 
n  n  n  n 

matrix  entries  M. .  and  the  entries  V.  for  i,j  <3  are  exactly  those  of 
in  i 

J  n 

(2-11).  To  obtain  the  other  M  _  ,  we  can  avoid  a  lot  of  writing  by  observing 
the  following:  each  row  is  augmented  by  the  coefficients  of  d^  and  e^, 
which  are  the  same  as  the  corresponding  coefficients  of  a^  and  b^  except 


n 

n 

replace 

j  by  n 
n  n 

in  Mn 
il 

to 

obtain 

(i  =  1, 

—  ,  6)  ; 

replace 

1  by  n 
n  n 

in 

i2 

to 

obtain 

■$5 

(i  =  1, 

. . . ,  6)  . 

Also,  equations  4,  5,  and  6  are  the  same  as  1,  2,  3,  except  that  they  are 
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For  example 


M45  =  (n2+n)  a(y2  nn  (y2)  “  nn  (y2))  ' 

A  number  of  degenerate  cases  can  be  obtained  from  these  equations  by  let¬ 
ting  the  appropriate  p  or  c  go  to  00  or  0 ,  and  thereby  eliminating  one  or 
more  unknowns  and  equations. 

Degenerate  '  !  ! 


Condition 

Shell 

1  Core 

'  Equations  (i) 

1  Unknowns  (i) 

- 

Elastic 

Fluid 

1 ,2 , 3 ,4, 5,6 

1 

1,2, 3, 4, 5, 6 

P4  ^  °o  ! 

Elastic 

l  Hard 

1 

1,2, 3/5,6 

'  1,2, 3,4, 5 

o 

II 

o. 

Elastic 

Soft 

1,2, 3, 4, 6 

1,2, 3, 4, 5 

o 

it 

X 

Elastic 

No 

1,2,3 

1,2,3 

o 

n 

CN 

u 

Fluid 

Fluid 

1,2, 4, 5 

1,3, 4, 6 

1 

c2=0,  P^00 

Fluid 

Hard 

1,2,5 

1,3,4 

c2=0,  P4=0 

Fluid 

Soft 

1,2,4 

1,3,4 

c2=0,  b=0 

Fluid 

l 

No 

1,2 

1,3 

P1  ■*  00  '• 

Hard 

No 

2 

3 

px  o 

Soft 

:  No 

1 

1  m-m  -  . 

3 

It  is  simplest  to  solve  this  set  of  equations  by  Gauss  elimina¬ 
tion  and  back  substitution.  Care  must  be  taken,  however,  since  the 
(complex)  entries  in  Mn  vary  in  magnitude  from  10^  to  10  Successful 

solution  is  achieved  by  normalizing  M  both  by  rows  and  by  columns  before 
applying  the  elimination  procedure.  The  same  considerations  hold  here  as 
discussed  in  Section  2.1  with  regard  to  the  order  of  the  Bessel  functions. 


Here  the  problem  is  somewhat  more  severe  because  of  the  larger  range  of 
Bessel  function  arguments.  In  the  examples  to  follow,  responses  will  be 
computed  to  ka  =  20  and  smoothed  impulse  responses  will  be  computed  for 
pulsewidth  parameter  a^  =  4. 


5.2  THIN  SHELL  LIMIT 


The  solution  for  the  thin  shell  can  be  obtained  nicely  from  the 
above  without  further  modification.  As  the  thickness  h  =  a-b  becomes 
small,  greater  accuracy  can  be  obtained  by  noting  that 


n  n  ,  3  n 

M  .  =  M.  .  +  h  -r—  M  . 
1]  4]  3r  43 


(j  =  1,  2,  4,  5)  . 


Row  4  can  then  be  substracted  from  row  1.  The  same  can  be  done  with  rows 
2  and  3.  As  h  +  0,  some  elements  of  Mn  can  be  replaced  by  zero.  However, 
this  is  not  a  great  simplification  since  no  unknowns  are  eliminated.  In 
the  calculations  that  follow,  the  unmodified  thick  shell  computer  program 
was  used  for  the  thin  shell  as  well. 


5.3  SHELL  RESPONSES  AS  FUNCTIONS  OF  ELASTIC  PARAMETERS 

The  sequence  of  plots  in  Figure  5-1  through  5-3  gives  the  mag¬ 
nitude  of  the  frequency  response  of  the  spherical  aluminum  shell  with 
parameter  thickness.  The  thickness,  h/a  decreases  from  1  to  .001  through 
the  sequence.  The  sequence  5-4  through  5-6  gives  the  same  results  in  the 

time  domain  for  smoothed  impulse  incidence  with  a  =4.  It  is  seen  that 

n 

the  thick  shell  has  a  more  complex  response  than  the  solid.  However,  as 
the  shell  becomes  thin,  the  response  becomes  simpler  again,  reducing  in 
the  limit  to  a  sequence  of  spikes  in  the  frequency  domain.  The  magnitude 
of  the  response  is  proportional  to  h  in  the  thin  shell  limit. 


Whereas  the  response  of  a  solid  is  perhaps  best  understood  in 
the  time  domain  as  a  sequence  of  pulses  corresponding  to  various  paths. 


H(ka)  H(ka)  H(ka) 


FIG.  5-1  Shell  responses  as  a  function  of  thickness  -  frequency  domain, 
thick  shell. 


the  very  thin  shell  response  is  best  understood  in  the  frequency  domain 
as  a  resonance  with  its  overtones .  The  resonant  frequency  is  dependent 
on  both  c^  and  c^  as  illustrated  in  Figure  5-7,  where  one  of  the  elastic 
parameters  of  an  aluminum  shell  with  h  =  .001  is  changed  at  a  time.  If 
both  and  h  are  changed  in  such  a  way  that  the  surface  density  (p^h) 
remains  constant,  then  the  response  remains  about  constant  (for  the  very 
thin  shell) .  This  is  illustrated  in  Figure  5-8. 

In  all  of  the  above  examples,  the  core  was  water  filled.  It 
is  interesting  to  see  the  effect  of  a  different  core.  In  Figure  5-9  the 
response  of  a  soft  sphere  is  compared  to  that  of  a  soft  core  (that  is,  a 
gas-filled)  spherical  shell.  The  difference  between  the  two  responses  is 
the  same  sequence  of  resonances  already  found  for  the  water  filled  shell. 

In  addition,  a  large  very  low  frequency  resonance  is  observed.  The  time 
domain  response  is  given  in  Figure  5-10.  A  similar  result  is  found  if  the 
core  is  a  heavier  fluid.  Figure  5-11  illustrates  the  response  of  a  glycerin 
core,  which  has  a  density  about  1.25  times  water.  This  resembles  closely 
the  response  of  a  fluid  sphere,  but  with  the  spikes  characteristic  of  the 
shell  superposed.  The  very  low  frequency  spike  is  absent. 

It  is  significant  that  a  thin  shell,  surrounding  a  sound-soft 
or  fluid  target  has  a  large  effect  on  the  response  in  the  form  of  sharp 
resonances . 
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(a)  SOFT 


FIG.  5-10  Comparison  of  responses  of  soft  sphere  and  soft  core  (gas-filled) 
spherical  shell  -  time  domain. 
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CORE  p  -  1.25 
FREQUENCY 
DOMAIN 


(b)  CORE  p-  1.25 
TIME  DOMAIN 


FIG.  5-1 1  Frequency  and  time  domain  responses  of  heavy  core  (glycerin-filled) 
spherical  shell. 
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SECTION  7 
APPENDICES 


7.1  ELASTIC  CONSTANTS 

Collected  here  are  some  useful  relationships  and  equations  to  do 
with  elastic  properties. 

X,y  Lame  Coefficients 
y,B,U  Young's,  Bulk,  Shear  Moduli 
0  Poisson  Ratio 

C1,C2  longitudinal ,  Transverse  Speeds 
p  Density 

Y  =  3B(1  -  0)  *  2y(l  +  o) 

X  =  B  -  2  y/3 
2  X  +  2y 

C1=“"F~ 


2  m  y 


Cl/C2  =  2(1  "  °)/(1  "  20) 


In  a  fluid: 


0,  0  =  *5,  <=2  =  0 


u  displacement 
Equation  of  motion: 

u  »  c*  7(7  •  u)  -  V  x  7  x  u 
In  fluid:  p  excess  pressure, 


p  .  -  c*  p(7  •  u)  ,  Vp  -  -  P  u 
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7.2  SCATTERING  MODEL  -  COMPUTER  PROGRAM 
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